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ABSTRACT 

This thesis discusses a procedure that optimizes the signal-to- 
noise ratio (SNR) detected by a linear array ina directional ambient 
noise field. The SNR can be optimized by minimizing the ambient 
noise detected by the array. For a given target location, each 
possible heading of the array centers the ambiguous beam of the 
array at a different true bearing. Therefore, each heading of the 
array will receive a different ambient noise level. An optimal 
heading can be obtained which maximizes the SNR received by the 
array. For all possible headings the beam pattern of the array must 
be determined and combined with the ambient noise field to 
determine the noise level detected. This thesis discusses the theory 
involved in calculating the ambient noise levels detected for each 


heading and provides a computer program which performs the 


calculations for a particular array. 





Thesis 
Mi ISH 


The reader is cautioned that computer programs developed in this 
research may not have been exercised for all cases of interest. While 
every effort has been made, within the time available, to ensure that 
the programs are free of computational and logic errors, they cannot 


be considered validated. Any application of these programs without 


additional verification is at the risk of the user. 
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I. INTRODUCTION 


The collapse of the Soviet Union and the distribution of diesel 
submarines to many third world countries have significantly changed 
the Antisubmarine Warfare (ASW) problem for the United States and 
its allies. The open ocean, deep water ASW scenario versus Soviet 
nuclear powered submarines has been replaced by what most 
experts agree to be a shallow water scenario against quiet diesel 
powered submarines. Additionally, improvements in submarine 
sound silencing throughout the world have drastically reduced 
passive acoustic detection ranges. The change in the ASW problem 
and the improvements in sound silencing on submarines has forced 
the United States to investigate new methods to increase passive 
detection ranges while maintaining the tactical advantage. 

The success or failure of a tactical towed array to detect a 
target of interest often depends on the signal-to-noise ratio (SNR) 
received at the sensors making up the array. Any opportunity to 
increase the SNR should be taken advantage of to increase the 
probability of detection. One viable means of increasing the SNR is to 
reduce the ambient noise detected by the array. The ambient noise 
detected by the array can often be reduced by exploiting the 
ambient noise field directionality by orienting the array in the 


optimal direction. For a suspected target location, each array heading 


centers the ambiguous beam (See Section IVA) at a different true 





bearing. Therefore, each heading will receive a different ambient 
noise level. By combining the beam pattern of the array with the 
ambient noise field for all possible array headings, an optimal 
heading can be obtained which minimizes the noise level detected by 
the array. | 
The majority of environmental models existing in the fleet 
today assume an omnidirectional ambient noise field. In many areas 
of the ocean the omnidirectional model is sufficient, while in other 
areas, including those close to shore and shipping lanes, the ambient 
noise level varies significantly with true bearing in the frequency 
spectrum dominated by distant shipping. The directionality of 
ambient noise can be used to the tactical advantage of the operator 
by orienting the linear array in a manner which aligns the 
ambiguous beam so that the total ambient noise received on beams 
directed towards the suspected location of the target is minimized. 
Linear array beam pattern calculations and the interaction of these 
patterns with the ambient noise field are discussed in detail in this 
thesis. Additionally, a computer program is provided which 
performs all the necessary calculations for a particular array. The 
program is limited toa single frequency and the target's suspected 
location must be known prior to running the program. The program 


is provided to demonstrating the feasibility of constructing a 


program for each tactical towed array in the Navy's arsenal. 














Il. HORIZONTAL DIRECTIONALITY OF AMBIENT NOISE 


A. SOURCES OF DIRECTIONAL AMBIENT NOISE 

The horizontal directionality of ambient. noise is caused by 
numerous sources including rain storms, biologics, and distant 
shipping. Most of these sources have short propagation ranges and 
rarely affect the directionality of ambient noise detected by passive 
acoustic sensors except for distant shipping. The frequency range in 
which distant shipping is the principle source, and therefore subject 
to high levels of directionality, is from about 50 to 500 Hz. [Ref. 1]. 
Within this frequency range the ambient noise can be dominated by 
noise which can travel for distances of 1000 miles or more from the 
location of the sensor by way of refracted raypaths in the deep 
sound channel. In this case, the ambient noise often has a distinct 
directional character with noise levels varying between true 
bearings. It is also within this frequency range that many passive 
acoustic sensors can detect enemy submarine signatures including 
blade rate information. 

There are two major mechanisms which account for the sound 
of distant shipping being trapped in the deep sound channel and 
therefore contributing to the directional character of the ambient 
noise field. The first mechanism is the tendency of the deep sound 


channel to approach the surface as the latitude increases [Ref. 2]. As 


can be seen in Figure 1, the depth of the deep sound channel in the 





Pacific Ocean ascends toward the surface as the latitude approaches 
the poles. This phenomenon occurs in the Northern Atlantic and 
Pacific where shipping traffic tends to concentrate when following 
great circle paths across the oceans. The large amount of shipping 
combined with the shallow deep sound channel in these areas cause 
a directional noise field in the vicinity of a passive acoustic array ~ 
which is located in the sound channel. This phenomenon accounts for 

the higher levels of directional ambient noise fields in the higher 


latitudes. 





Figure 1. Depth of SOFAR axis in the Pacific Ocean. Contour interval 
200 meters [Ref. 3]. 


The second mechanism accounting for shipping noise being trapped 


in the deep sound channel is the "downslope conversion 








phenomenon" [Ref. 2]. This describes the process of bottom 
reflected/surface reflected raypaths changing to refracted raypaths 
as the sound crosses over the continental shelf into the deep ocean 
basin and into the deep sound channel. Dense shipping normally 
found close to shore causes this phenomenon to cause significant 
. directional noise towards the shore when operating off the 
continental shelf. Figure 2 displays both of these mechanisms 
occurring at the location marked on the chart. The lobe to the 
northwest can be attributed to the high levels of shipping traffic 
located south of the Aleutian Islands traveling along the great circle 
track between the Far East and North America. The lobe towards the 
northeast may be caused by the down-slope conversion of the sound 
of shipping on the continental shelf off the west coast of the United 
States [Ref. 4]. As the ASW mission moves from an open ocean, deep 
water scenario to a scenario closer in to shore, the down-slope 
conversion phenomenon will become a more important factor 
contributing to the directional nature of the ambient noise field. 
B. OBTAINING THE DIRECTIONAL NOISE FIELD 
The directionality of the noise field can be obtained from noise 
modeling based on historical data or from real time in-situ 
- measurements. The in-situ measurements have the highest levels of 
confidence of being valid for a given scenario. 
1. Noise Modeling 


The amount of noise generated by distant shipping is 


directly proportional to the amount of shipping traffic present. | 






















= 


Figure 2. Horizontal Directionality of Ambient Noise. Site location 
marked by X on top. Directional patterns, as obtained with end-fire 
beam of line array (dashed) and with broadside beams (solid) on 
bottom [Ref. 2]. 
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Therefore, the directionality of the ambient noise field can be 
estimated by the amount and location of the shipping in the vicinity 
of the sensor. Historical data are available on the shipping density 
around the world and can be used to model the expected directional 
noise field for a given location at a given time of year. The 
environmental factors affecting sound transmission in the oceans 
must be modeled as well. 
2. Noise Directionality Measurements 

The most reliable way to obtain the directionality of the 
ambient noise field is to measure it with on board sensors. This can 
be done by performing any number of different towship maneuvers 
that will allow the measurement of beam noise data on three 
different array headings. In order to resolve the ambiguities, each 
heading must be at least forty degrees apart. Wagstaff's Iterative 
Technique [Ref. 5] is an ambiguity resolution and noise-field 


estimation algorithm which can be used to perform in-situ 


calculations of the directional noise field. 











Ill. TOWED ARRAY BEAM PATTERN 


A. BEAMWIDTH 
All equations discussed below are taken from Lawrence J. 
Ziomek's book Underwater Acoustics [Ref. 6]. 
1. Beamwidth at Arbitrary Steering Angle 
When a given towed array beam pattern is steered from 
broadside towards end-fire, the beamwidth increases. A minimum 
beamwidth value occurs at the broadside beam providing the finest 
bearing resolution. To compute the beamwidth, first assume that the 
array is linear and is oriented along the X axis. Next, let us calculate 
the beamwidth for the horizontal beam pattern. As seen in Figure 3, 
w' is the steering angle in degrees, measured counter clockwise from 
the positive X axis. Also, let ¥+ and W- represent the angles 
corresponding to the 3-dB down points of the horizontal beam 
pattern. Note that Figure 3 shows both of the half-beamwidth angles 
to be equal to AW/2. Although this is not true in general, the 
resulting equation for AW will still give accurate results. 
In the XY plane the direction cosine u with respect to the 
X axis is given by 
u = cos (W), (1) 


and will be used for calculation purposes. Therefore, let 
u+ = cos(¥+) >O (2) 
u- = cos(W-) >O_ (3) 











Up>u, (4) 
Wy = UW - (AW/2), (5) 

and | 
w= w+ (AW/2), (6) 





x 


Figure 3. Representation of the beam steering angle W' and the 
beamwidth AY [Ref. 6]. 


Thus, 
Au = Uy - u- = 2sin(#")sin(AW/2) > O (7) 
since 
cos(a-B) - cos(a+p) = Z2sin(a)sin(p) (8) 
_ giving 
AW = 2sin-1[Au/2sin(¥")], (9) 


where AW is the beamwidth in degrees of the horizontal beam 
pattern centered at the angle ¥'. Note that wp = n/2 corresponds to 
the broadside beam pattern, and that W' = O or x corresponds to the 
end-fire beam pattern where Au/2sin(W') = © and Eq. (9) is not valid. 
Equation (9) is only valid when lAu/2sin(W')| <1. Therefore, a 


different approach must be used to determine the beamwidth at 


end-fire geometry. 











2.  Beamwidth at End-Fire Geometry 
Restricting ourselves to the horizontal beam pattern, 


consider the case when the steering angle, W', is zero. Letting 


u- = cos(¥-) > O (10) 
and substituting Y'=0 into Eq. (6) yields 
W. = AW/2. (11) 
For this case, 
u' = cos (#") = cos (0) = 1 (12) 
and 
u- = 1 - (Au/2). (13) 
Combining Eq. (13) with Eqs. (3) and (11), . 
cos(AW/2) = 1 - (Au/2) (14) 
or 
AW = 2cos-l(1 - (Au/2)) (15) 


where AY is the beamwidth in degrees of the endfire beams. 
B. BEAMPATTERN 

Consider a linear array composed of N (odd) identical, equally 
spaced, amplitude weighted elements. For the purposes of this 
thesis, only rectangular amplitude weights will be discussed. It is 
important to note that the theory pertains to all types of amplitude 
weighting. However, by limiting the discussion to rectangular 


amplitude weights, a closed form expression can be obtained for the 


directivity function of the array. 








1. Directivity Function 
The directivity function of a linear array can be 
expressed as a function of the frequency and the appropriate 
direction cosine, which in our case is direction cosine u. Figure 4 
shows a typical normalized beam pattern steered in the direction u = 


u' in direction cosine space. 


ID (fa)! 





Figure 4. Magnitude of a typical normalized beam pattern steered in 
the direction u = u' in direction cosine space [Ref. 6]. 


The closed form expression for the normalized, horizontal directivity 
function of a linear array with rectangular amplitude weighting is 
Dy (ffx) = sin(xfxNd)/Nsin(xfxd) ' (16) 
where 
fx = u/A, (17) 





N is the number of elements in the array, d is the inter-element 
spacing, u is given by Eq. (1), and Ais the wavelength of the sound 
waves. Upon steering the beam, the resulting directivity function, 
Dy'(f,fx), is simply a translated version of the original directivity 
function, Dy(f,fx), that is, 
Dn'‘( ffx) = Du(ffx-fx’). (18) 
Therefore, when steering the beam, Eq. (16) can still be used by 
replacing fx with fx-fx' where: 
fy-fx' = (u-u')/2. (19) 
he Beamwidth in Direction Cosine Space 
To determine the 3-dB beamwidth, Au, in direction cosine 


space for the normalized directivity function given by Eq. (16), set 


sin(xNdu,/a)/Nsin(xdu,/d) = 0.707 (20) 
or, rewriting, | 
sin(Nzx) - N(O.707)sin(«x) = 0 (21) 
where 
x= u,d/A= Aud/2), (22) 
so that | 
Au = 2Ax/d. (23) 


Upon finding the solution to Eq. (21) for x and substituting the result 


into Eq. (23), the 3-dB beamwidth in direction cosine space is 
obtained. 





TV. OPTIMAL LINEAR ARRAY HEADING 


The towed linear array is one of the most valuable sensors on 
board any Naval ASW unit. The high levels of directivity obtained 
and the relatively high degree of mobility combine to make the 
towed linear array an effective ASW sensor. These two factors have 
also led to the common practice of placing the narrow broadside 
beam in the suspected direction of the target. This tactic has worked 
well in the fleet, especially when the ambient noise field is 
omnidirectional. However, in many cases where the ambient noise 
field is directional, this tactic can lead to a severe degradation of the 
towed array's performance. The interaction of the ambiguous beam 
with the ambient noise field can lead to a less than optimal noise 
level received by the array. 

A. AMBIGUOUS BEAMS 

The beamwidth of a linear array beam pattern has a minimum 
value at broadside and increases to a maximum value at the endfire 
beams. This conical beam pattern inherently presents a left-right 
ambiguity problem as seen in Figure 5. It is impossible to determine 
from which side of the array the sound is coming without 
maneuvering the array. Each beam of the array has a mirror image 
on the opposite side of the array commonly called the ambiguous 
beam. The ambiguous beam receives additional noise which, in 


many cases, raises the total noise level received by the array which 
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centers the main lobe of the beam pattern on the target. It is often 
possible to minimize the ambient noise arriving at the ambiguous 


beam by prudently choosing the heading that the array is oriented. 


Aft 
"Endfire” 
Beam 


Array 
Centerline 


weeeece 





Figure 5. Generic Linear Array Beam Pattern. 


B. TOTAL NOISE LEVEL CALCULATION 
When numerous individual noise sources are presented to a 
linear array, they must all be considered in determining the total 
noise level. Assuming that the individual noise levels are 
uncorrelated, the total noise level may be calculated using the 
following equation: | 
TNL = 10 log(3,10NLi/10)) (24) 





where NLj are the individual noise levels in dB and TNL is the total 
noise level [Ref. 7]. 
C. OPTIMAL HEADING DETERMINATION 

Figure 6 displays a directional ambient noise field in dB ref 1 
uPa in 1 degree increments. The target's suspected location is to the 


northwest. 


Target 










Ambiguous Beam 


Main Beams for for Heading 1 


Headings 1 and 
2 


Ambiguous Beam 
for Heading 2 


Array a 


Heading 2 _ 
Array 
Heading 1 


Figure 6. Illustration of the interaction of array beams and the 
ambient noise field. 


Array heading 1 is to the south and places the target in beams about 
45 degrees off of broadside (solid beam pattern). The ambiguous 


beam on heading 1 is located to the northeast where the ambient 
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noise averages about 50 dB. Using Eq. (24) and neglecting detailed 
beam pattern calculations, the approximate total noise level received 
for heading 1 is about 62 dB. Array heading 2 is to the southwest 
and places the target in the broadside beam (dashed beam pattern). 
The ambiguous beam for heading 2 is located to the southeast where 
the ambient noise level averages about 60 dB. Again using Eq. (24) 
and neglecting detailed beam pattern calculations, the total noise 
level received for heading 2 is about 67 dB. Array heading 1 
receives a total noise level which is about 5 dB less than that of array 
heading 2. The SNR received on array heading 1 will subsequently 
be greater than the SNR received on the broadside array heading 2 
and the probability of detection will likewise be proportionately 
higher. Assuming that the target's signal strength does not vary 
‘significantly over the time it would take to maneuver the array, 
heading 1 should be chosen over the broadside heading 2. 

This technique can be used to determine the actual noise level 
detected by the array, including detailed calculations for the beam 
pattern, for all possible headings. The total noise levels detected for — 
each heading can be compared with the noise levels for all the other 
possible headings to determine the heading which detects the 


minimal amount of noise. Thus, an optimal heading can be 


determined for the orientation of the array which maximizes the SNR 
detected by the array [Ref. 8]. 





V. PROGRAM DESCRIPTION 


The program in appendix A is provided to demonstrate the 
manner in which the optimal heading is obtained for a given 
directional noise field and for a suspected target location. 

A. AMBIENT NOISE DATA 

The program prompts the user to input the ambient noise field 
by one of two means. First, the user can input the name of a file 
where the ambient noise data is stored. This would correspond toa 
~ historical type of ambient noise data which has been stored in the 
computer at a previous time. The second means of inputting the 
ambient noise data is for the user to input the data directly into the 
computer. This process is somewhat cumbersome and takes about 
five to ten minutes depending on the complexity of the noise field, 
but the data is copied to a file for easy access at a later time. This 
corresponds to an in-situ ambient noise field measurement. 

B. LINEAR ARRAY CHARACTERISTICS 

The program simulates a linear array comprised of eleven 
omnidirectional hydrophones with an interelement spacing of five 
meters. The interelement spacing allows for optimal array 
performance for 150 Hz signals, which is the only frequency the 
program is capable of processing. The number of elements generate 


a broadside beam which is nine degrees wide and an endfire beam 


which is forty-six degrees wide. The characteristics of this array do 





not correspond to any particular array, but are used purely to 
demonstrate that any array can be simulated with a similar program. 
C. TARGET INFORMATION 

The only information that is required concerning the target is 
the target's suspected true bearing. The user is prompted to input 
this information which the program uses to determine the relative 
bearings to the target for all possible headings that the array may be 
oriented. Additionally, the program determines the headings which 
place the broadside beam in the direction of the target. The 
broadside headings are used as a reference later in the program. 
D. BEAM PATTERN DETERMINATION 

For each possible heading, the main lobe is centered at the 
suspected location of the target. The program also determines the 
bearing at which the ambiguous beam is centered for each heading. 
The 3-dB beamwidths of the beams are determined using either Eq. 
(9) or Eq. (15), depending on the orientation of the target to the 
array. The program also determines the normalized beam pattern 
using Eq. (16). The true bearings which the main lobe and the 
ambiguous beam cover are determined for use in the noise 
calculation. 
E. NOISE LEVEL CALCULATION 

The program combines the beam pattern and the ambient noise 
field to determine the ambient noise level detected for each possible 
heading. The ambient noise level for each bearing which the beam 


covers is converted into a pressure, multiplied by the beam pattern 
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factor, then summed together with the other noise levels in the same 
beam pattern using Eq. (24). After all the pressures are summed the 
total noise level is converted back into a dB level for later use. For 
each heading, the total noise level detected is compared to the - 
reference broadside heading total noise level. The results are stored 
in a file which can later be displayed using another program. 
F. PROGRAM LIMITATIONS 

The program is capable of simulating one linear array and can 
only process noise data at 150 Hz. Additionally, only rectangular 
windowing is used when determining the beam pattern of the array. 
Future possible work could provide programs which can simulate 
multiple arrays and frequencies. Other types of windowing could 
also be incorporated into the program. 
G. PROGRAM EXAMPLE 

Figure 7 displays an ambient noise field at 150 Hz given in dB 
referenced to 1 wPa. The noise levels are given in one degree 
increments in degrees true. The higher noise levels are concentrated 
to the south while lower noise levels are found to the north. If the 
suspected location of the target is to the north, the broadside 
headings of 090 or 270 would place the main broadside lobe in the 
direction of the target. While these courses place the main lobe in 
the direction of the target, the ambiguous beam is directed towards 
the south in the direction of the higher ambient noise levels. 
Therefore, courses which still place the main lobe in the direction of 


the target, but place the ambiguous beam in areas of lower noise 
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levels may have lower total noise levels than the broadside headings. 
Even though the beamwidth of the array beam pattern is greater 
than the broadside beamwidth and more individual noise levels will 
be detected, the total noise level may be less than the broadside 
heading total noise level. 

The optimal heading program was run for the case described 
above and the results are displayed in Figure 8. The SNR gain for 
each possible heading relative ‘to the broadside heading are 
displayed. The possible headings are given in degrees true and the 
SNR gains are given in dB. Note that for the broadside headings of 
090 and 270, the SNR gain is zero relative to broadside which is 
expected. As the array heading is changed from the broadside 
heading towards the north, the total noise level received by the 
array lowers and the SNR gain relative to the broadside heading 
increases. Note that when the array is oriented about 30 to 60 
degrees off of broadside in any direction, the SNR gain is from 5 to 6 | 
dB relative to the broadside heading. For this case the common 
practice of placing the broadside heading in the direction of the 
target degrades the performance of the array by about 6 GB. 
Additionally, the heading to orient the array on can be chosen to 
comply with any other constraints placed on the ASW unit. 

Additional examples using various ambient noise fields are 


presented in Appendix B along with a discussion of each one. 


21 











to the broadside 


headings relative 


Figure 8. SNR gain for all true 


heading. 





VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The optimal array heading can, in many cases, significantly 
improve the SNR for a linear array in a directional noise field. Two 
major applications for the technique occur at different phases of a 
given ASW scenario. The first application occurs when planning to 
search an area for a target. The traditional barrier search is a prime 
example of this. The target's suspected location is known and there 
is plenty of time to consider the directional ambient noise field prior 
to setting up the barrier. The orientation of the barrier should be 
based upon the optimal array heading calculation for the given noise 
field. 

The second major application occurs after a target has been 
detected. The ASW unit should place itself between the target and 
the direction of the highest ambient noise to prevent 
counterdetection from the target submarine while optimizing the 
unit's detection capabilities. Additionally, the ASW unit should 
choose the courses which optimize the array's SNR based on the 
optimal heading calculation. . 

B. RECOMMENDATIONS 
Based on the results of this thesis, the following actions are 


recommended: 
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1. A computer program, such as in Appendix A, should be 
made up for each of the arrays in the US Navy's arsenal and 
distributed to the applicable ASW units capable of using them. 

-2. A fleetwide ambient noise directionality measurement 
procedure should be made up and distributed to the fleet. All ASW 
unit commanders should be advised to the practicality of measuring 
the ambient noise field and using the results to optimize the 
probability of detection. | 


3. A directional ambient noise data base should be established 


which can be made available to all ASW units. 
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APPENDIX A 
FORTRAN PROGRAMS 


This appendix provides the reader with the program "Optimal" 
which is compiled using the FORTRAN language. The program 
calculates the optimal heading to search for a target in a directional 
ambient noise field. The program is described in detail in Chapter V. 
Additionally, two plotting programs are provided for display 
purposes. The first program, "Plot", displays the directional ambient 
noise field using the Disspla routine on the Naval Postgraduate 
School's mainframe computer. The second program, "Plot2", displays 


the results of the Optimal program using the Disspla routine also. 
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FILE: OPTIMAL FORTRAN Al 


DAVID C. MCDONNELL, LT, USN 


3 HHA IRIE TEI IEIEIDE HE HIER IE IESE SEE IIE HE HEHE TIE HE IEEE IE IE HEE NE TE IE HE FEE HEI 
XTHIS PROGRAM CALCULATES THE OPTIMAL TOWED ARRAY HEADING IN X 
¥A DIRECTIONAL NOISE FIELD ENVIRONMENT, MINIMIZING THE % 


XAMBIENT NOISE DETECTED BY THE SENSOR. % 
HEH HEHE HEHEHE HE DEE DE IK HE HE HE IE DE DE 33 HE HE HE 3 3 3 3 EE HE IE HE 3 3 DEE HE HE HE BE 3E DE IE 3E HE HE HE DE DE DE IE 3K 3 DE 3 HE HE IE IEEE HE 


VARIABLE DECLARATION SECTION 


CHARACTERX20 FILE1,FILE2, FILES 

REAL X(C360),Y(360),TGT, (180), N1,N2, DN, HDG, DU,N, BH, DTHETA, DT1 
: REAL DT2, DT3, DT4, PSUM, PI,ANS, LEVEL, BHEAD, PSYIP, PSY2P, F,F1 

REAL G(180), PSY1, ae poate? DPSY2 

INTEGER A,B,C,D, eS »0,P,Q 


3 HEE HE IE HE IE HE HE DK HE HE IE HE 3 HE BE HE IE HE HE HE HE DE HE DE HE HE DE HE DE HE IE HE IE BE IE DE HE BE IE 3K DE HE 3 HE 3 SE IE 3E 36 DE HE IE IE 3E 3E 36 3E HE IEE 


VARIABLE DESCRIPTION SECTION 


x FILE1, FILE2, AND FILE3 ARE USER FILES TO ACCESS AND STORE 
INFORMATION. 

*X0360) IS AN ARRAY WHICH STORES THE AMBIENT NOISE DATA 
*¥C360) IS A DUMMY ARRAY. HOLDING THE AMBIENT NOISE DATA 

*TGT IS THE TARGETS TRUE BEARING 

¥LC€180)} IS AN ARRAY WHICH STORES THE AMBIENT NOISE LEVEL FOR 
EACH HEADING BETWEEN 1 AND 180 TRUE 

*¥N1 IS THE NOISE LEVEL AT THE BEARING OF THE TARGET 

*N2 IS. THE NOISE LEVEL AT’ THE AMBIGUOUS BEARING OF THE TARGET 
*DN IS THE DIFFERENCE IN BEARING BETWEEN TARGET AND SENSOR » 
*HDG IS THE HEADING OF THE SENSOR 

*DU IS A CONSTANT WHICH AFFECTS THE BEAMNIDTH OF THE ARRAY 

*N IS THE NUMBER OF ELEMENTS IN THE ARRAY 

XBW IS THE BEAMHIDFH OF THE ARRAY IN RADIANS 

¥DTHETA IS THE. HALF-BEAMWIDTH OF THE ARRAY IN DEGREES 
*DT1,DT2,DT3,DT4 ARE THE TRUE BEARINGS WITHIN THE MAIN BEAM 
*PSUM IS THE SUM OF THE NOISE PRESSURES 

xPI=3.141593 

¥ANS = ANSHER TO QUESTION 

XLEVEL IS THE NOISE LEVEL IN DB OF THE AMBIENT NOISE REF 1UPA 
XBHEAD IS THE BROADSIDE HEADING TO THE TARGET 

XPSY1P AND PSY2P ARE THE ANGLES IN DEGREES CCW FROM THE SENSOR 
HEADING TO THE TARGET AND THE AMBIGUOUS BEAMS 

*F IS THE BEAM PATTERN FACTOR 

¥Fl IS A_ FACTOR USED IN CALCULATING THE BEAM PATTERN FACTOR 
¥GC180) IS THE ARRAY WHICH STORES THE GAIN FOR EACH HEADING 
XPSY1,PSY2,DPSY1,DPSY2 ARE FACTORS USED TO CALCULATE THE BEAM 
PATTERN FACTORS 

¥A,B,C,D,1I,J,M ARE INTEGERS USED FOR NOISE FIELD INPUT 
¥K,0,P,Q ARE INTEGERS USED FOR THE GAIN CALCULATION 


PT=3.141593 
HH HEHEHE HE HEHEHE HE HEH HE HEHE HE HE HEH HE HEHE HE HEHE HEH HE HE HE HEHE HE BEE HE DE 3 HEHE DE HE HE HE DE 3 HE HE IE 3E IE 3E IE HE HE 36 3 IC HE IE TEE 


OBTAIN THE AMBIENT NOISE DATA FROM THE USER 


THE USER IS GIVEN A CHOICE TO USE EITHER A STORED DATA FIELD - 
OR TO CREATE HIS ONN DATA FIELD 


WRITECX,*) "PLEASE CHOOSE THE NUMBER CORRESPONDING TO THE * 
WRITEC*,*)'MEANS OF ENTERING THE AMBIENT NOISE FIELD* ) 
WRITECK,*)*  # 

HRITECX, x)! 1) OBTAIN DATA FROM AN EXISTING FILE! 

WRITECX, *)! 2) USER INPUT DATA 

READ(C¥,*,ERR=5) ANS 


NQAAANQAANND 


UMAOOANOAN ONNAAANAAAANNANAANANAAANANANNANANAANNAN|ANDNNN 


c 
c IF THE USER ene TO USE A STORED DATA FIELD 
. IFCANS.EQ.1) THE 
6 WRITECX, X) 'ENTER THE NAME OF THE FILE WHERE THE NOISE?® 


WRITECX,*) "DATA IS STORED! 
READC(, *CA10)", ERR= 6) FILE1L 
OPENCUNIT=1, FILE=FILE1, STATUS='OLD') 
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FILE: OPTIMAL FORTRAN Al 


DO 20 1=1,360 
READC1,%) XCID 





20 , CONTINUE 
CLOSE (1) 
GO TO 70 
c 
c IF THE USER WANTS TO EREATE HIS OWN DATA FIELD « 
ELSE IF CANS.EQ.2) THE 
D=1 
c Mel 
21 WRITE (%,%) ‘ENTER THE FILE WHERE YOU WISH THE NOISE* ° 


WRITE (%,%) "LEVELS TO BE STORED! 
READ (%,'(A10)',ERR=21) FILE2 ; 
OPEN CUNIT=2, FILE=FILE2, STATUS="UNKNOWN® ) 


c 
25 WRITE (%,26) D 
26 FORMATC*ENTER THE NOISE LEVEL BEARING!,1X,13.3,1X,"DEG TRUE") 
WRITE (€%,%)*CENTER -99 TO END DATA ENTRY WHEN REACH BRG 360)! 
c READ (*,%,ERR=25) LEVEL 
IF (LEVEL.£EQ.-99) GO TO 60 
29 WRITE (%,*) "ENTER THE BEARINGS WHICH CORRESPOND TO THIS LEVEL* 
WRITE (€¥,30) D 
30 FORMAT C*STARTING WITH',1X,13.3,1X,"'DEG TRUE CEX. 001,010)*3 
READ (%,*,ERR=29) A,B 
C=B-Atl 
D=B+1 
c » 
c THIS DO LOOP WRITES THE NOISE LEVELS TO THE FILE 
c 
DO 50 J=1,C 
WRITE (2,40) LEVEL 
40 FORMAT (1X,F6.2) 
XCMI=LEVEL 
MelH+1 
590 CONTINUE 
c GO TO 25 . 
690 CLOSEC2) 
c 
70 END IF 
Cc BE HE HE HEHE HE HE HE HE HE HE HEHE HE HE HE HE HE HE HE HE HE DE HE HE DE DE HK FE IE HE HE HE HE DE HE BE IE DE IE IE DE IE HE DK DE DE DE 3G HE BE DE IE 5E HE HK DE HE HE IE IE HE 
c HEHE HEH HE HE HEHEHE HK HEH HK HE HE HEH HEH HEHE HEHEHE HK HE HE HE HEH HEH HEH HEHEHE HK HEH HEH HEHEHE HK HHH HEHE HE HHH HEH HHH HK 
c . 
c THIS PROGRAM USES A GENERIC LINEAR ARRAY COMPOSED OF 7 
c ELEMENTS SPACED 5 METERS APART, GIVING THE ARRAY A LENGTH OF 
c 50 METERS. THIS SPACING IS SET UP FOR A FREQUENCY OF 150 HZ. 
c FOR PURPOSES OF BEAMWIDTH DETERMINATION THE FOLLOWING EQUATION 
& HAS PREVIOUSLY SOLVED FOR X. 
€ SINCNXPEXX)-N%. 7071 XSINCPIXX)=0 
c RESULTING WITH X = .04041, AND DU=4%X GIVING: 
c 
=.1616 
C N=11 
c EG I DE DE HE 3E IE IK 3G HE 3 DE IE 3 DE IE DE BE DK DE 3G HE 3 3E DE DE DE IE 3E DE DE 3G HE DE DE HE DE IE HE DE DE 3K 3 3E DE DE IE DE BE IE DE HE HK DE 3E IE EE EE 
& OBTAIN THE TARGET INFORMATION FROM THE USER 
75 PIRITEC%,*) "ENTER THE SUSPECTED TRUE BRG OF THE TARGET! ~ 


READCX,¥, ERR=75) TGT 
IF (TGT.EQ.0) TGT = 360. 


DETERMINE THE BROADSIDE HEADING BETWEEN 1 AND 180 FOR THE 
GIVEN TARGET POSITION. 


Hann 


BHEAD = TGT - 90 
IF (BHEAD.LE.0) THEN 

BHEAD = TGT + 90 
ELSE IF (BHEAD.GT.180) THEN 
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BHEAD = BHEAD - 180 


END IF 
EK HEHEHE HEH HEHEHE HEHEHE HEHE I I HE HE HE HEHE HE DE IE IE HE IE HE DE IE 3 BE 3 DE HE DE DE DE DE 3K 3 3 3 HE 9 9 3K 3 3 DE E26 HEE 3E 36 3 HE HEE 


c 
Cc 
Cc OBTAIN NAME FOR FILE WHERE SNR GAIN DATA IS TO BE STORED 
Cc 
80 WRITE €%,%)' ENTER FILENAME FOR SNR GAIN DATA STORAGE? 
READ (%,'CA10)',ERR=80) FILES 
OPEN CUNIT=3, FILE=FILE3, STATUS='UNKNOWN! ) 
. c 
WRITE €3,%) TGT 
c 
Cc 3G HE HEHE HE HE DE DE HE HE DE HE IE HE HE HK HE HE DE IE HE HE BE DE HE DE DE IE DE DE DE HE HE DE DE HE DE IE DE DE IE DE HE IE 3E IE DE DE DE DE I HE IE HE IE IE HE DE IE IE HEE 
- Cc . 
c FOR EACH POSSIBLE HEADING THE LOOP CALCULATES THE GAIN 
c RELATIVE TO BROADSIDE THAT THE ARRAY ACHIEVES 
c . 
HDG = 1 
c 
DO 200 K = 1,180 
c 
NL=TGT 
c 
DO 90 Q=1,360 
¥€Q)=XCQ) 
2° CONTINUE 
DN=HDG - TGT 
c 
c THE FOLLOWING IF STATEMENT DETERMINES THE AMBIGUOUS BEAM FOR 
c THE CURRENT HEADING 
c . . 
IF (DN.GT.-360.AND.DN.LE.-180) THEN 
DN=DN+360 -~ 
N2=HDG+DN 
ELSE IF (DN.GT.-180.AND.DN.LT.180) THEN 
N2=HDG+DN 
ELSE IF (DN.GE.180.AND.DN.LE.360) THEN 
DN=360-DN 
N2=HDG-DN 
c END IF 
IF (N2.GT.360) THEN 
N2=N2 ~ 360 
ELSE IF (N2.LT.0) THEN 
NH2=N2 + 360 
c END IF 
IF (N1.EQ.0) N1=360 
c IF (N2.EQ.0) N2=360 
Cc HHH HK HHH HHH HHH HEH HH HEHEHE HEH HEHEHE HEHE HE HEHE HEHEHE HE HE HE HE HE DE HE HE HE HE HE HEE HE HEHE HEH HEHEHE HK HK 
c 
C DETERMINE THE BEARING CCN FROM THE ARRAY FOR THE BEAM 
c CONTAINING THE TARGET AND THE AMBIGUOUS BEAM. THESE RELATIVE 
c BEARINGS ARE USED IN CALCULATING THE BEAMNIDTH AND THE BEAM 
C PATTERN FACTOR USED BELOW 
IF CHDG.GE.TGT) THEN 
- PSY1P=HDG-TGT 
PSY2P=360-PSY1P 
ELSE IF CHDG.LT.TGT) THEN 
PSYIP=HDG-TGT+360. 
. PSY2P=360-PSY1P 
c END IF 
HEHE I HE HE HE HE HH HE HK HE HE HEHE HE DE HE HE HE HE HE HE HE IKE HE HE IE IE IE FE 3 HE HEHE HE DE FE DE HE HE HE HE HE DE HE HE EE HE HE HE DE HE 3E IE HE IE IE HEE 
c DETERMINE THE BEAM WIDTH OF THE BEAM CONTAINING THE TARGET AS 
c WELL AS THE AMBIGUOUS BEAM. THE 3 DB DOWN POINT IS USED TO 
€ DETERMINE THE BEAMWIDTH OF THE. BEAM PATTERN. 
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Cc 
c 


Qa00 


2 OH20H0N ANNANDA AAANNDN O20 


FOR ENDFIRE 


IF CPSY1P.EQ.0,.OR.PSYIP.EQ.180.) THEN 
BH = 2x%ACOS(1-(DU72)) 


ELSE IF NOT ENDFIRE 


ELSE IF CPSYIP.NE.0..AND.PSYIP.NE.180.) THEN 
IF CDUr COMABSCSINCPSYIPXPT7180))). LE.1) THEN 
BH = 2*ASINCDU/C2¥ABSCSINCPSY1PXPI7180)))) 
ELSE IF (DUs C2XABS(SINCPSY1P¥PI/180))).GT.1) THEN 
ean ite BH = 2*ACOSC1-(DU/2)) 


END IF ; 

DETERMINE THE HALF BEAMWIDTH OF THE BEAM IN DEGREES. 

DTHETA = INTCBWX1807(PIX2)) , 

HE HE HE IE HE HE HE 3E BE HE HE HE FE IE HE HE HE IE DE HE HE DE 3E HE 3E HE DE DE DE DE IE HE DE HE IE HE IE HE HE HE HE IE IE DE HE DE DE HE HE HE IE IE HE HE HE IE DE DE HE HE ETE 
CALCULATE THE TOTAL NOISE SEEN BY THE ARRAY ON THIS HEADING 
PSUM = 10X*CXCINTONI))710) + 10%%(XCINT(N2) 710) 

THE FOLLOWING DG LOOP SUMS THE AMBIENT NOISE LEVELS WITHIN THE 
BEAMS OF THE ARRAY. THE BEAM PATTERN FACTOR, F, IS CALCULATED 
FOR EACH NOISE LEVEL en THE BEAM. 

DO 100: 0=1,DTHETA 


NOISE LEVELS ARE SET TO ZERO AFTER COUNTING TO PREVENT THEM — 
FROM BEING COUNTED TWICE. 


DTl=Wl1+0 © 


IF (DT1.GT.360) DT1 = DT1 - 360. 
PSYL=PSY1PXPI7180. 
DPSY1=CPSY1P+0)*PI/180. 
F1l=CPI72)*CCOSCPSY1)-COSC¢ DPSY1)> 
IF ae -EQ.0) THEN 


LS 

F=SINCNXF1)/(NXSINCFL)) 
END IF 
IF CF.LT.0.707) GO TO 92 


PSUM=PSUM+F%10**CYCINTCDT1))710) 
YCINTCDT1)9=0 


DT2 = N2 + 0 
IF (DT2.GT.360) DT2 = DT2 - 360. 


PSY2=PSY2PXPI/180. 
DPSY2=(PSY2P+0)*PI/180. 
Fl=CPI72)%(COSCPSY2)-COSCDPSY2)) 
IF aks THEN 


LSE 
F=SINCN®F1)/(NXSINCF1)) 
END IF , 
IFCF.LT.0.707) GO TO 94 


PSUM=PSUM+F®10%*CYCINTCDT2))710) 
YCINTCDT2))=0 


DT3 = Nl 
IF CDT3. be: 0 DTS = DT3 + 360 


DPSY1=CPSY1P-0)XPI/180. 
FLl=CPI72)*(COS(PSY1)-COS(DPSY1)) 
IF CF1.EQ.0) THEN 
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F=1. 
LSE 

F=SINCNXF1)/ CNXSINCF1)) 
END IF 
IF CF.LT.0.707) GO TO 96 


Cc 
PSUM=PSUM+FX10%*CYCINTCDT3))/10) 
YCINT(DT3)9=0 
c 
. 96 T4 = N2 
IF (DT4.LE. 8) DT4 = DTG + 360 
Cc 
DPSY2=(PSY2P-0)*PI/180. 
- F1=(PI/2)*C(COS(PSY2)-COS( DPSY2)) 
IF (F1,EQ.0) THEN 
F=1. 
LSE 
F=SINCNXF1)/ CNXSINCFL)) 
END IF 
c IF (F.LT.0.707) GO TO 100 . 
PSUM=PSUM+F¥1 0X*CYCINTCDT4))/10) 
C YCINTCDT4)3=0 
100 CONTINUE 
116 L(K)=10*L0G10CPSUM) 
HDG = HDG +1. 
200 CONTINUE . 
c THE FOLLOWING DO LOOP CALCULATES THE GAIN RELATIVE TO , 
C BROADSIDE FOR EACH HEADING. 
C DO 250 P = 1,180 
c FOR EACH HEADING THE GAIN RELATIVE TO THE BROADSIDE BEAM IS 
¢ CALCULATED AND PLACED INTO THE ARRAY GCP). 
C GOP) = LCINTCBHEAD)) - LCP) 
WRITE (3,%) GCP) 
250 CONTINUE 
C THE FILE WHERE THE GAIN DATA IS STORED IS NOW CLOSED 
CLOSE (3) 
c 
END 
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c DAVID C. MCDONNELL, LT USN PL000010 
c THIS PROGRAM USES THE DISSPLA ROUTINE TO PLOT THE AMBIENT NOISE  PLOO0020 
c FIELD BEING USED. THE USER MUST IDENTIFY THE FILE WHERE THE NOISEPLO00030 
c DATA IS STORED. A POLAR PLOT IS USED TO DISPLAY THE NOISE LEVEL Pt0o00040 
c IN DEGREES TRUE. : PLO00050 
c : PLO00060 
c VARIABLE DECLARATION SECTION PL000070 
c PLOOO080 
INTEGER I,J PLO00090 ~ 
REAL X(360),Y¥(360),Z2(360),MIN, MAX, DR, DRI, DR2, DRS PLOOO100 
CHARACTERX20 FILES,FLNS PLO0G110 
c PL000120 
c OBTAIN FILE NAME FROM THE USER Froagis9 - 
WRITE (%,%) ENTER FILENAME AND FILETYPE WHERE DATA IS STORED'* PLOO0150 
READ (%,'C€A20)}*)} FILES _ PLGO00160 
c PLO00170 
c A "7" MUST BE APPENDED TO THE FILENAME FOR THE DISSPLA ROUTINE PLO001380 
c TO RECOGNIZE THE FILE. PL000190 
c PLOO0200 
FLNS="/°//FILES//°At PLO00210 
c PLO00220 
- OPEN CUNIT=1,FILE=FLN$,STATUS="OLD!) BL ogacse 
c THE FOLLOHING DO LOOP READS THE AMBIENT NOISE FROM THE FILE. THE PLOOGO250 
c DISSPLA COMMAND "POLAR" PLOTS THE DATA REFERENCED TO RADIANS FROM PLO000260 
c EAST. THEREFORE, THE DEGREES TRUE ARE TRANSFORMED INTO RADIANS ‘PLO00270 
c FROM EAST FOR PLOTTING PURPOSES. PLO00280 
c PLO00290 
c INITIALIZE VARIABLES PL000300 
MIN=100. : . PLODO3S10- = 
MAX=0. PLO00320 
c PL000330 
BO 10 I=1,360 oe PLO00340 
READ (1,%) XC€I) PLO00350 
c : PL000360 
c DETERMINE THE MINUMUM AND MAXIMUM NOISE LEVEL FOR PLOT PURPOSES ae 
IF (XCI).LT.MIN) MIN=XC€I) ; PLO00390 
IF (XCI).GT.MAX) MAX=XCI) PLO00400 
c PLO00410 
IF CI.GT.O.AND.I.LT.90) THEN PLG00420 
J=90- PL000430 
ELSE IF (CI.GE.90) THEN PL000440 
| J=450-I PLO00450 
END IF PLO00460 
YCJ=REALCJ) PLO00470 
Z(5)=XCI) PLO00480 
c PLO00490 
10 CONTINUE PLO00500 
Cc PLO000510 
c SET MIN/MAX EQUAL TO +7-10 THAN THE ACTUAL MIN/MAX NOISE LEVEL PLO00520 
c FOR PLOTTING PURPOSES. PLO000530 
MIN=INTCMIN)-20. PLO00540 
. MAX=INTCMAX)410. PLO00550 
| DR=INTCCMAX-MIN) 7/6.) PLO00560 
DRI=MIN+DR PLO00570 
DR2=DR1+DR PLO00580 
DR3=DR2+DR PL000590 
DR4=DR3+DR PLO06600 - 
DR5=DR4+DR PLOO0610 
MAX=DR5+DR PLO000620 
c . PLO00630 
DR=2*DR PLOG0640 > 
c PLOO00650 
c PLO00660 
c THE PLOTTING ROUTINE USES DISSPLA COMMANDS TO PLOT THE DATA PLOO0670 
c IN DEGREES TRUE WITH THE AMBIENT NOISE LEVELS GIVEN IN DB REF PLO00680 
c 1 UPA. PLO00690 
c PL000700 
CLOSECUNIT=1) PL000710 
CALL SHERPAC'*MCT3NOIS', 'A',3) PL000720 
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FILE: 


c 


PLOT 


FORTRAN Al 


CX41 (4107) 
PAGE (8.5,11.) 
AREA2D(6.,6.) 

HEADIN (* AMBIENT NOISE 
POLARC .01745,DR,3.,3. 
POLORG (MIN) 
GRID(-30,2) 
MARKER (15) 
CURVECY,Z,360,1) 
REALNOCMIN, +0,2.8 
REALNOCDR1,+0,3.3 
REALNOCDR2,+0,3.8 
REALNOCDR3,+0,4.3 
REALNOCDR4,+0,4.8 
REALNOCDR5,+0,5.3 
REALNOCMAX,+0,5.8 
MESSAGC *NOISE LEV 
MESSAGC'N',1,3.0, 
MESSAGC'180',3 
MESSAGC 090° ;3 
MESSAGC '270',3 
ENDPL (0) 
DONEPL 


NENTS SS 


GON Oevuvuvvrn 


7 


PLOT’,20,1.,1) 


VEN IN DB REF 1UPA* 
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734,1.0, 


~1.0) 


PLO000730 
PLO00740 
PLO00750 
PLO00760 
PLO00770 
PLOO00780 
PL000790 
PLO00800 
PLOO0810 
PLO00820 
PLO000830 
PLO00840 
PLOO0850 
PLO00860 
PLO00870 
PLOO0880 
PLO00890 
PLOO0900 
PLOOO0910 
PLOO9920 
PLO000930 
PLO00940 
PLOO0950 
PLOO0960 
PLO00970 
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MAANAHANAHD 


an00 


amAag O00 


ielelele) 


ann 


aaeMnn 


o 


AQAQAN FO 
o 


DAVID C. MCDONNELL, LT USN 

THIS PROGRAM USES THE DISSPLA ROUTINE TO PLOT THE AG DATA 
GENERATED FROM THE OPTIMAL PROGRAM. THE PROGRAM READS THE 
DATA FROM THE FILE AND PLOTS THE DATA _ON A POLAR PLOT. — 
THE SHIPS HEADINGS AREREPRESENTED IN DEGREES TRUE AND THE 
AMPLITUDE OF THE AG IS REPRESENTED IN DB REF 1 UPA. 


VARIABLE DECLARATION SECTION 

INTEGER I,J,K 

REAL X(360),Y(360),2(360),A,B,MIN, MAX, DR, DR1, DR2, DR3, TGT 
REAL DR4,DR5,STEP,T1 

CHARACTER®20 FILES, FLNS 

OBTAIN THE FILE NAME FROM THE USER WHERE THE DATA IS STORED 


WRITEC%, ¥) "ENTER FILENAME AND FILETYPE WHERE DATA IS STORED* 
READ(*,*(A20)') FILES 


A ¢ MUST BE APPENDED TO THE FRONT OF THE FILENAME FOR THE 
DISSPLA PROGRAM TO RECOGNIZE IT. 


FLNS="/'//FILES//°At 
INITIALIZE VARIABLES 


A=0. 
MIN=100. 
MAX=0. 


THE FOLLOWING DO LOOPS READ THE DATA FROM THE FILE. THIS IS DONE 


THICE DUE TO THE RECIPROCAL HEADINGS WILL HAVE THE SAME AG. 


DO.20 I=1,2 ~ om 
OPEN CUNIT=1,FILE=FLN$,STATUS="OLD') 
READ(C1,*) TGT 
IF (TGT.LE.90) T1=90-TGT 
IF (TGT.GT.90) T1=450-TGT 
DO 10 J=1,180 
A=At+1] 
READ (1,%) XCINTCA)) 


DETERMINE THE MINIMUM AND MAXIMUM VALUE OF DATA FOR PLOT PURPOSES 


IF CXCINTCA)).LT.MIN) MIN=XCINTCA)) 
IF CXCINTCA)).GT.MAX) MAX=XCINTCA)) 


THE DISSPLA POLAR COMMAND PLOTS ANGLES IN RADIANS FROM EAST. 
THEREFORE DEGREES TRUE ARE CONVERTED TO RADIANS FROM EAST IN 
THE FOLLOWING IF STATEMENTS. 


IF Saenger Eeagne THEN 
ELSE IF €A.GE.90) THEN 
B=450-A 
END IF 
YCINTCB) )=REALCB) 
ZCINTCB)SXCINTCAD) 


CONTINUE 
CLOSE (1) 
CONTINUE 


SET MIN/MAX EQUAL TO +7- 10 FOR PLOTTING PURPOSES 


MIN = MIN - 10. 

MAX = MAX + 10. 

DR=INTC CMAX~MIND 76 3 
STEP=2xDR 

IF Sede nOE tT SeeR)? THEN 


DR3=DR4-DR 
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PLOGO0010 
PL000020 
PL000030 
PLO00040 
PLOCO050 
PLO00060 
PL000070 
PLOO0080 
PLO00090 
PLO00100 
PL000110 
PLO00120 
PLO00130 
PLO00140 
PLO00150 
PL000160 
PLO00170 
PLO00180 
PLO00190 
PLO00200 
PLO00210 
PLO00220 
PLO00230 
PLO00240 
PLOOO250 
PLO000260 
PLO00270 
PLOOG280. 
PLOOG290 
PLO00300 
PLOG0310 
PLO000320 
PLO00330 
PLO00340 
PLO00350 
PLO00360 
PLO00370 
PLO00380 
PLO00390 
PL000400 
PL000410 
PLO00420 
PLQG0430 
PLO00440 
PLOGO45D 
PLO00460 
PLO004790 
PLO000480 
PLO00490 
PLO00500 
PLOOGG510 
PLO00520 
PLO00530 
PL000540 


.PLO000550 


PLO00560 
PLO00570 
PLO00580 
PLG00590 
PLO00600 
PL000610 
PLOO0620 
PLO00630 
PLO00640 
PLO00650 
PLO00660 
PL000670 
PLO00680 
PLO00690 
PL060700 
PLOOO0710 
PL000720 


FILE: PLOT2 FORTRAN Al 


DR2=DR3-DR PL000730 
DR1=DR2-DR PLO00740 
MIN=DR1-DR PLO00750 
DR5=DR4+DR PLO00760 
MAX=DR5+DR PLO00770 
ELSE IF (CO-MIN).LT.(C3%DR)). THEN PL000780 
DR2=0 PLO00790 
DR1=DR2-DR PLOO0800 
MIN=DR1-DR PL0008190 
_ DR3=DR2+DR PLO00820 
DR4=DR3+DR PLO00830 
DR5=DR4+DR PLO00840 
MAX=DR5+DR _ PLOO0850 
- END IF . PL000860 
Cc . PLG00870 
c THE PLOTTING ROUTINE USES DISSPLA COMMANDS TO PLOT THE DATA PLO00880 
c IN DEGREES TRUE WITH LEVELS GIVEN IN DB REF 1UPA. Matt ttt 
c : : 

CALL SHERPAC *MCTSDATA', tA*,3) PLo00910 
c CALL CX41(4107) PL000920 
CALL PAGE (8.5,11.) PL000930 
CALL AREA2DC6.,6.) PLO00940 
CALL HEADIN C*SNR GAIN RELATIVE TO BROADSIDE HEADING’, 38,1.,1) PLO00950 
CALL POLARC.01745,STEP,3.,3.) PLO00960 
CALL POLORG CMIN) PLO000970 
CALL GRID (-30,2) PLO00980 
CALL MARKER(C4) PLO000990 
CALL CURVEC(T1,MAX,1,-1) PL001000 
CALL MARKER (15) PL001010 
CALL CURVE CY,Z,360,1) PL001020 

CALL REALNOCMIN,+0,2.8,2.7) PLO01030 » 
CALL REALNOCDR1,+0,3.3,2.7). PLO01040 
CALL REALNOCDR2,+0,3.8,2.7) PLOO01050 
CALL REALNOCDR3,+0,°4.3,2.7) . PLO01060 
CALL REALNOCDRG,+0,4.8,2.7) PLOO01070 
CALL REALNOCDR5,40,5.3,2.7) PLO61080 
CALL REALNOCMAX,+0,5.8,2.7) PLO01090 
CALL MESSAGC 'LEVELS GIVEN IN DB REF 1UPA*,28,1.3,-1.0) PLOO1100 
CALL MESSAGC'N',1,3.0,6.1) PLOO1110 
CALL MESSAG(C'180',3,2.8,-.2) PLO01126 
CALL MESSAG('090!,3,6.2,3.0) PLOO1130 
CALL MESSAGC'270',3,-.5,3.0) PL001140 
CALL MESSAGC('X = TARGET TRUE BEARING’, 23,1.3,~-1.4) PLOO1150 
CALL REALNOCTGT,+0,4.2,-1.4) PLO01160 
CALL ENDPL (0) PL001170 
CALL DONEPL PLOO1180 
SToP PLOO1190 
END PLO01200 
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APPENDIX B 
ADDITIONAL EXAMPLES 


Examples of the "Optimal" program, along with a discussion of 
each, are included in this appendix. The first three examples are test 
cases which partially validate the program, and the other examples 
show the applicability of the program, depicting the SNR gain 
achievable in some cases when optimal heading selection is used. 

A. EXAMPLE 1 

Example 1 uses an omnidirectional ambient noise field which 
has a level of 50 dB as seen in Figure 9. In this case the broadside 
heading should be the optimal heading due to the narrow beamwidth 
of the broadside beams. Figure 10 shows the results of the "Optimal" 
program for this case with the suspected target location bearing 090. 
Notice the highest SNR gain is for the headings oriented towards the 
north and the south, the broadside headings. As the linear array's 
' orientation is changed away from the broadside headings, the SNR 
gain becomes negative. Thus, the array's performance is degraded by 
up to about 4 dB when the array is oriented to place the target in the 
endfire beams. When the array is oriented to place the target just 
off of endfire, the array's performance is even worse than that of the 
endfire headings. This is due to the fact that the main lobe and the 


ambiguous beam overlap and a large heart shaped beam pattern is 


formed which has a larger beamwidth than that of the endfire beam. 




















B. Example 2 

Example 2 provides another test case which is presented to 
partially validate the "Optimal" program. An ambient noise field 
which has a very high noise level at just one true bearing, 090, was 
created and is displayed in Figure 11. With the suspected target 
location to the west, the broadside headings of 000 and 180 would 
place the ambiguous beam in the direction of the high noise level at 
090 and the array performance would be severely degraded. Any 
heading which does not place the array's ambiguous beam in the 
direction of the high noise level will have a significant SNR gain 
relative to the broadside heading. As can be seen in Figure 12, the 
SNR gain for almost all headings is anywhere from 11 to 18 dB 
higher than the SNR for the broadside heading. Figure 12 shows the 
results of the array's performance as the ambiguous beam passes 
through the high noise level located at 090. Note the curved pattern 
of the SNR gain plot in Figure 12 showing the inverted beam pattern 
of the array for the headings which place the ambiguous beam in the 
high noise level area. Additionally, total noise level calculations 
using Eq. (24) verify the program's result of the 18 dB SNR gain for 
the headings just off broadside. | 
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SNR GAIN RELATIVE TO BROADSIDE HEADING 


| 





LEVELS GIVEN IN OB REF IPA 
X = TARGET TRUE BEARING 270. 


Figure 12. SNR gain for all true headings relative to the broadside 
heading in a highly directional noise field. Note the inverted beam 
pattern for array orientation of 000 and 180. . 
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C. EXAMPLE 3 . 
Example 3 is another highly directional ambient noise field 
which is used to partially validate the "Optimal" program. In this 
case very high noise levels of 80 dB, compared with the rest of the 
noise field, are located between true bearing of 160 and 200 as seen 
in Figure 13. With the suspected target location towards the north, 
the broadside headings of the array place the ambiguous beam in the 
direction of the high noise levels. The array's performance will be 
degraded for headings which place the ambiguous beam in the high 
noise level area. Figure 14 shows the results of the "Optimal" 
program for this case. The SNR gain for all headings which do not 
place the ambiguous beam in the high noise level area perform much 
better than the broadside headings of 090 and 270. Again the 
inverted beam pattern of the array is represented centered about 
the broadside headings. The constant SNR level of 0 dB gain results 
when the entire ambiguous beam is in the high noise level area. As 
_ the ambiguous beam starts to exit out of the high noise level area, 
the SNR gain gradually increases up to 27 GB, at which point none of 
the ambiguous beam is in the high noise level area. Total noise level 
calculations using Eq. (24) verify the program's result of a 27 GB gain 
for the headings just off broadside, which do not place the ambiguous 


beam in the high noise level area. 


42 


AMBIENT NOISE PLOT 








NOISE LEVELS GIVEN IN DB REF 1JPA 


Figure 13. Highly directional ambient noise field at 150 Hz. Levels 
are given in dB ref 1 uPa. Note the high noise level located towards 











Figure 14. SNR gain for all true headings 
heading in a highly directional noise field. 
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D. EXAMPLE 4 

Figure 15 displays an actual ambient noise field. Noise levels 
vary from a high of about 60 dB towards the east-southeast to a low 
of about 45 dB towards the west. The "Optimal" program was used to 
analyze for a suspected target location towards bearing 315. As seen 
in Figure 16, gains of up to 3 dB can be obtained by choosing the 


appropriate array headings in this case. For instance, choosing a 


array heading of about 285 will produce a SNR gain of 3 dB relative 
to the broadside heading of 045. 
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E. EXAMPLE 5 

Figure 17 is an example of an actual ambient noise field which 
has two high noise level peaks of 60 dB at true bearings of about 050 
and 085. Additionally, the noise levels toward the west-southwest 
are lower with levels of only about 45 dB. The "Optimal" program 
was run for this case and the results are displayed in Figure 18. Note | 
the large variation in the SNR gain with relatively small changes in 
array heading. For instance, for an array heading of 330, a 4 dB SNR 
gain is achieved, but for an array heading of about 335, the SNR gain 
is -2 CB. A 6 dB difference in array performance is obtained with 
only a change of 5 degrees in array orientation. The two peaks in the 
ambient noise field are responsible for the high variability in the SNR 


gain in this case. 
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SNR GAIN RELATIVE TO BROADSIDE HEADING 








LEVELS GIVEN IN DB REF i1pPA 
_ X = TARGET TRUE BEARING 270. 


Figure 18. SNR gain for all true headings relative to the broadside 
heading in an actual directional noise field. 
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